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Abstract 

We calculate order a s color magnetic corrections to the valence 
quark distributions of the proton using the Los Alamos Model Poten- 
tial wavefunctions. The spin-spin interaction breaks the model SU(4) 
symmetry, providing a natural mechanism for the difference between 
the up and down distributions. For a value of a s sufficient to produce 
the N — A mass splitting, we find up and down quark distributions in 
reasonable agreement with experiment. 

1 Introduction 

Historically, quark models have provided us with a convenient, if simplistic, 
method for making quantitative calculations of low energy hadronic observ- 
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ables. Unfortunately, much of the available data on hadronic structure is 
taken in an energy regime well beyond that at which simple quark models 
are expected to be valid. In particular, the distribution of quarks in the nu- 
cleon as a function of their light cone momentum fraction has been measured 
in numerous experiments over a wide energy range in the scaling regime. In 
order to take advantage of this data, it was necessary to invent an argument 
allowing an extrapolation from the low energy quark model regime to the 
high energies where direct measurements can be performed. 

Such an argument was put forward by Jaffe and Ross|L| in 1980. Quark 
models, they argued, could only be a representation of QCD at a relatively 
low renormalization scale, /i 2 ~ ^/Rp ro ton> as a resu lt of their relative simplic- 
ity. More explicitly, they argued that, at large /x 2 , experiments are able to 
resolve structures of order and consequently one sees a hadron as a very 
complicated object, composed of valence quarks, lots of glue, and a sizeable 
sea. As the renormalization scale gets smaller, the resolution decreases, sea 
quarks get recombined into gluons, and gluons are reabsorbed by quarks. 
Indeed, one can speculate that at very low renormalization scales, virtually 
all of the sea and most of the glue can be reabsorbed into the valence quarks, 
so that one is left with a proton composed of three "constituent" quarks 
bound together by some residual interaction, a quark modelpj! While this 
argument says nothing about which of the great variety of quark models de- 
scribes the nucleon, it provides a well-defined prescription for testing models 
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against high energy data. One simply calculates quark distribution functions 
using the model and then uses the renormalization group to evolve the result 
from the quark model scale, /Xq, to a scale relevant to experiment. 

The first step is to calculate the quark distribution at the quark model 
scale. For unpolarized scattering, the spin averaged quark distributions can 
be written, in the nucleon rest frame, as[||] 

<&(*) = ^ / dtS^t- < NMCh + P + (C)MO)\N > \lc 

(1) 

where q + = —Mx/^/2, with x the Bjorken scaling variable, iftifyi) are quark 
field operators of flavor i, j + = (7 + 7 3 )/v^2 is a Dirac gamma matrix, 
P±(£~) is a path ordered exponential, exp[±ig s Jq A + (r]~)dr]~] with g s the 
strong coupling constant, that insures the gauge invariance of the operator, 
and the subscript LC denotes the light cone condition on £, namely £j_ = 
£ + = 0. The approach we shall take has been described in detail in Refs. 
3-5 for several models, and consists of a straightforward evaluation of the 
matrix elements in Eq. 1, using a Peierls-Yoccoz momentum eigenstate|| to 
describe the nucleon in its rest frame. While other approaches exist in the 
literature || , our belief is that the current method maintains closer contact 
with the quark model in question, and is more easily generalized to systems 
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of more than one nucleon, which we will examine elsewhere. 

Phenomenologically, a problem arises when the quark distributions are 
evaluated using the unperturbed, SU(4) symmetric nucleon wavefunction. 
As a result of the symmetry, the u and d valence quark distributions have 
the same functional form, and consequently d v (x)/u v (x) = 1/2 for all x. 
Experimentally, it is well known that the ratio decreases as x increases, and 
is thought to vanish linearly as x — > 1. The missing ingredient, it seems, is a 
mechanism for SU(4) breaking. 

Fortunately, two mechanisms immediately present themselves. The first 
is the idea that the nucleon is surrounded by a cloud of virtual mesons, mainly 
pions, and that the pions contribute a flavor asymmetry to both the valence 



and sea distributions |10[. While this idea is certainly worthy of study, it is 
not the subject of the current effort. Instead, we shall concentrate on the 
second mechanism for SU(4) breaking, the color magnetic interaction. 

The essential idea here is that the spin-spin interaction results in a depen- 
dence of the quark wavefunctions on the spin state of the other quarks that 
reside in the nucleon. Since, in the naive SU(4) wavefunction of the nucleon, 
the spin is unevenly distributed among the different flavor quarks, the spin 
dependence of the color magnetic interaction is transmuted into a flavor de- 
pendence of the spin averaged quark wavefunctions. Put into the calculation 
of quark distributions, this flavor dependence in the wavefunctions turns into 
a flavor dependence of the quark momentum distributions. 
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In this paper, we calculate the color magnetic corrections to the valence 
distributions of the nucleon using the Los Alamos Model Potential(LAMP). 
In the next section, we briefly describe the LAMP, and calculate its quark 
distributions in the absence of color magnetic effects. We then begin the 
calculation of the color magnetic corrections, separating out contributions to 
the valence distributions involving two and three body effects, respectively. 
We then describe the contributions of the gauge correction to the valence 
distributions, and how they may be partially resummed to obtain improved 
distributions. Finally, we shall evolve the full expressions for the corrected u 
and d distributions to high Q 2 , where they will be compared with experiment, 
and discuss the results. 

2 Two and Three body Corrections 

In order to set the stage for the forthcoming 0(a s ) calculations, we begin 
by calculating the unperturbed quark distributions for the LAMP model]!]]]. 
The model consists of three massless quarks, bound together in a linear scalar 
potential, 

V{r) = V {r-r ), (2) 

with parameters Vq= 0.9 GeV/fm and ro=0.57 fm. Wavefunctions and single 
particle energies are obtained by solving the Dirac equation for this poten- 
tial. As is usual for quark models, the potential parameters are chosen to 
generate the average of the nucleon and delta masses, and the color magnetic 
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interaction is invoked to generate the nucleon-delta splitting. 

In Ref. 5, the matrix elements in Eq. 1 are evaluated, assuming that the 
nucleon is described by a Peierls-Yoccoz momentum eigenstate, 

\N,P = 6) = \J d 3 a\N, R cm = a), (3) 

where A is a constant required to covariantly normalize the state, and 

— * — * 

\NR cm = a) denotes the unprojected state with center fixed at R cm = a. If 
one assumes that the time dependence of the quark field operators is well 
approximated by the the single particle eigenvalue of the Dirac equation, the 
expression for the valence quark distribution is given by 



M = -^-[[J \ dkG{k){tl{k)+t\{k) + -^t Q {k)t 1 {k))] 



2k_ 

dkU[k)[tn[k) +tUk) + 

+[*-->*+]!, (4) 



where 



G(k) = Jrdr sin kr A 2 (r)EB(r) , 

V = J r 2 drA 3 (r)EB(r), 

k(k) = Jr^rUkrMr), 

ti(k) = J r 2 drji(kr)v(r), 

A(r) = Jd 3 z4(z-^M^, 



(5) 



with ipo(r) the ground state, single quark wavefunction, with upper and lower 
components u(r) and ia ■ fv(r)/r, k± — uj ± Mx, with u the single particle 
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energy, and EB(r) = (EB,R cm = r\EB, R cm = 0) is the matrix element 
between two "empty bags" separated by a distance r, which accounts for 
the dynamics of the confining forces. In the LAMP picture, the dynam- 
ics of the confining fields are not specified, so we have no guidance on the 
proper choice for this function||. For the remainder of this paper, we assume 
EB(r) ~ const. The resulting distribution is plotted, along with those of 
the MIT[|J and soliton J7] bag models in figure 1. The LAMP distribution 
lies between the MIT bag, where the valence quarks carry all of the mo- 
mentum at the bag scale, and the soliton bag, where roughly a quarter of 
the proton's momentum is carried by the confining degrees of freedom. We 
expect, therefore, that the quark model scale, liq, from which we need to 
evolve the LAMP distribution in order to compare with experiment, will lie 
somewhat below the .6 GeV 2 found in reference 5 for the soliton bag. If we 
vary our assumption of a constant EB(r), it is easy to generate distributions 
that interpolate continuously between the curve shown in figure 1 and curves 
resembling the soliton bag distribution. Correspondingly, the scale /i 2 , will 
increase to a value in the neighborhood of .6 GeV 2 . 

Before calculating the color magnetic corrections to Eq. 4, we must en- 
gage in a brief digression on the questions of renormalization prescriptions 
and gauge choices. As it stands, Eq. 4 represents the results of a QCD in- 
spired model of the structure of the nucleon. If the model were, in fact, a 
solution of the QCD Hamiltonian in some calculational scheme, the choice of 
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gauge and renormalization scheme would be specified, and the quark model 
wavefunctions would be unambiguously defined. Unfortunately, there is no 
known scheme for generating the LAMP (or any other) model directly from 
QCD, and consequently the gauge and renormalization schemes, if any, that 
produce the wavefunctions of reference 9 are unspecified. In order to com- 
pletely specify the model, it is necessary to postulate that the model wave- 
functions are calculated in a particular gauge, with a particular choice of 
renormalization scheme. As these choices are necessarily ad hoc, it follows 
that different choices will, in general, result in different models, and neces- 
sarily different predictions for physical observables. For the purposes of this 
paper, we shall choose the MS renormalization scheme, and A = gauge. 

The strategy we shall employ for evaluating the color magnetic correc- 
tions to the valence quark distributions is a straightforward application of 
perturbation theory, augmented by a closure approximation whose scale is 
the mean energy of the gluon in the nucleon. The interaction Hamiltonian 
is given by 

H I = g 9 J d 3 rA a (r,t)J a (r,t), (6) 

where g s is the strong coupling constant, A a (f,t) is the gauge field, and 
J a (r,t) = ^jj(f,t)^^ip(f,t) is the colored quark current operator, with 4r an 
SU(3) generator, and 7 a Dirac matrix. To leading order in g s , the nucleon 
wavef unction is given by 

I AT, p = 0) = Z-*[\N , P = 0) + 1 Hj\N , P = 0)], (7) 

hj — n 
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where Z is a wavefunction renormalization constant, and E is the ground 
state energy before correction. The energy shift is given by 

6E = ^V(iVo, P = 6\H IT ^-—H I \N , P = 0), (8) 

with \4o the volume of space, and the renormalization constant Z is 

Inserting a complete set of gluon states, and a complete set of colored, 3 
quark states, and using the fact that the unperturbed Hamiltonian from Eq. 
2 is independent of color, the energy denominators may be seen to be given 
by 

H -E = E R {k) + u{k) - Eo, (10) 
where k is the gluon momentum, u(k) its energy, and En(k) is the energy 

— * 

of an excited 3 quark state R with momentum —k. In order to parallel the 
assumptions regarding the gluon propagator in reference 9, and to facilitate 
evaluation of the numerical integrations to come, we make the ansatz that 
u(k) = /ie fc2 / 2At2 , where \x is an effective gluon mass taken to be 400 MeV. 
The energy denominator is then given by 

H - E = /ie fc W(l + e-*'*(E R (%) - E )/fi). (11) 

A reasonable expectation is that the OGE process couples dominantly to 
the ground state, and that the relevant gluon momenta are small (on the 
order of 1/ R p ) , so that the recoil energy of the massive nucleon is small, and 
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the second term in parentheses can be neglected. With this assumption, the 
energy denominator becomes independent of the intermediate three quark 
excited state, and closure may be invoked to perform the sum over excited 
states. The results of equations 7-9 keep the same form , except that the 

— * 

factors E — H are replaced by the constant u(k). 

Setting P ± = 1, the 0(a s ) corrections to eq. 1 may be obtained using 
eq. 7. The result is 

Sq(x) = -Z- 1 ^ J dCe' tMx ^^ J d 3 rd 3 r'G(f-f) 

x (N , P = 0| J a (r>(f)7 + ^(0) J o (0 \N , P = 0) \ LC 

+ {Z~ l - l)q°(x), (12) 

where J a (r) = 4>(r)?fip(r), G(r) = Ae~^ r2 l A is the effective gluon propa- 
gator in the closure approximation, with A oc a s = |^ an overall strength 
parameter chosen to reproduce the nucleon-delta mass splitting, and q°(x) is 
the unperturbed quark distribution of eq. 4. 

2.1 Two-Body Correction 

In general, eq. 12 gives rise to three types of terms which may be classified 
by the number of valence quarks acted on by the three currents appearing in 
the matrix element. The "one-body" terms, in which all three currents act 
on the same quark, represent self interactions and are traditionally ignored 
in quark model calculations. Here, these corrections will tend to mimic the 
effects of QCD evolution, and may be absorbed into a redefinition of the 
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scale where the evolution begins. Since a precise determination of this 
parameter is not the purpose of the present effort, we shall honor tradition 
and neglect the one-body terms. 

Next, there are the "two-body" terms, which arise when a gluon is in the 
process of being exchanged between two valence quarks when one of these 
quarks is struck by an external probe, such as a photon. Finally, there are 
three-body terms, in which the two spectator quarks exchange a gluon while 
the struck quark interacts with the external probe. These terms cannot be 
absorbed into a redefinition of /x , since as we shall see, they are dependent 
on the isospin of the struck quark. 

Using fermion anticommutation relations, the two and three- body correc- 
tions to the valence quark distributions can be extracted in a straightforward 
manner. The result for the two-body correction is 



where a denotes the struck quark, (3 the spectator with whom the gluon is 

— * 

exchanged, ipo(k) is the quark momentum space wavefunction, the brackets 
denote an average over the color wavefunction of the proton, and 




x 



(13) 
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with 

Zp{ri,r 2 ) = J AG(f)^(£-fi)7^(f-f 2 ). (15) 

In order to extract the physically relevant u and d quark distributions, we 
must average over the spins in eqs. 13-15. To do the average, we must separate 
the spin dependence of the quark wavefunctions appearing in eq. 13. For 
s-wave quark wavefunctions, we get 

4>oa(ki)^°7 + il>oa(k 2 ) = hC Q1 (k 2 , h) - k 2 C Q1 (k 1 , k 2 ) 

+i((k 2 x <7 Q )C i(£;i,£;2) + (h x a a )C 01 (k 2 , h)) 
+z(C 00 (k 1 , k 2 ) - (ki ■ k 2 + ik 2 ■ (h x a a ))Cu{ki, k 2 )) 
+iz x ((T Q (C o(/ci, k 2 ) + h ■ k 2 C u {ki, k 2 )) 
-(k 2 fa ■ a a + k ± k 2 ■ <j a + i(ki x fc 2 ))Cn(A;i, fc 2 ))),(16) 

where 

^^-ttW 1 - (17) 

and G (0, 1) as defined in eq. 5. The function Zp(ki, k 2 ) can be decom- 
posed similarly, 

Zp{k u k 2 ) = k x Z X p{k u k 2 ) + k 2 Z 2 p (k u k 2 ) 

x 3 p )Zl{k 1 , k 2 ) + (k 2 x a^Z^h, k 2 )), (18) 

where the Zp(ki, k 2 ) are scalar functions of k\ and k 2 , evaluated semi-analytically 
using eq. 15 and the six gaussian fit to the quark wavefunctions of reference 
9. 
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Once the integrations are performed, there are only three vectors left, 
z,a a and ap, and the two-body correction takes the form 

p+o. 

where T 2b (x) is a tensor defined in terms of the quark wavef unctions. 

Using rotational arguments, it is easy to demonstrate that only the trace 
of T will contribute to the spin independent quark distribution. Hence, the 
functions f 2b (x) and Tu(x) may be written as 

f ( x ) = T~i — / /* \r, o(Mx-u-z- fci) 



3/i 2 J (2tt) 5 
x Z 2 (fci, k 2 ) (kjdoih, k 2 ) + h ■ hCoiih, k 2 ) 

+z ■ ki(C 00 (ki, k 2 ) + ki ■ k 2 C u (ki, k 2 )fj 

+Z 1 (k 1 ,k 2 )(k 1 ■ k 2 C 10 (h,k 2 ) + klCQi(k u k 2 ) 

h ■ k 2 



+z ■ ^(J—^Cooih, k 2 ) + klCuih, k 2 ))) 



(20) 



and, 



—5(Mx -w-z-ki) 



(2tt) 



x 



Z 4 (k 1 ,k 2 )( K k 2 1 C 10 (k 1 ,k 2 ) - k x ■ k 2 C i(k 1: k 2 ) 
+z ■ fci(Coo(A;i, k 2 ) - ki ■ k 2 Cn(k u k 2 fj 
+Z\k 1 ,k 2 ){k l -k 2 C w {k 1 ,k 2 )-k 2 2 C m {k 1 ,k 2 ) 
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+z ■ k 2 (C o(h, k 2 ) - h ■ k 2 Ci 1 (k 1 , fc 2 )) , 

(21) 

Assuming that the unperturbed proton wavefunction is SU(4) symmetric, 
the spin sums are given by 

E *« • = "3 ± 1, (22) 

so that the two-body corrections to the physical quark distributions are given 
by 

8u 2b (x) = 2f 2b (x)-2T 2b (x) (23) 
Sd 2b (x) = f 2b (x)-AT 2b (x). (24) 

An interesting feature of this result is that the d quark distribution is twice 
as sensitive to color magnetic effects than the u distribution. On physical 
grounds, this is understandable, since the two spectator u quarks are neces- 
sarily in a spin 1 state in the SU(4) wavefunction, while the spectator ud pair 
is most likely to be found in a spin state. Numerically, the color electric 
correction, f 2b (x), is quite small and is neglected in the following. 

The correction to the antiquark distribution is obtained from eq. 13 by 
replacing u <-> — u and including an overall minus sign. 

2.2 Three Body Correction 

After separating out the one and two-body terms from equation 12, only the 
three-body term remains. Physically, this correction accounts for the mod- 
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ified structure of the spectator diquark, and is expressed as a modification 
of the recoil function G(k) appearing in equation 4. In terms of the quark 
wavefunctions the three body correction is given by 



x$ 0a (ty-f + 'ip 0a (k)Fp e (-k), 



d 3 k Mx — uj — z • k 
( 2?r M 7k ) 



V2 



(25) 



where the function Fp e (k) is given by 



F Pe (k) = / d z y^F p€ {y)EB{y), 



(26) 



with 



Fp e (z) = J J dVidV 2 (7(ri-f 2 )Vty3(n)7V^^^ (27) 

Just like the two-body correction, Fp e {z) can be decomposed into a scalar 
function, f 3b (z) and the trace of a tensor T 3b (z), 

(28) 



F Pe (z)=f\z) + ^T 3b (z), 



where 



f 3b (z) = -2 J J d 3 r 1 d 3 r 2 G(f 1 - r 2 ) v(r 1 )v(r 2 )u(f 1 - z)u(f 2 - z)f\ ■ f 2 



-v(ri)v(f 2 - z)u(r 2 )v(f 2 - z)f 1 



(r 2 - z) 



\r 2 - z\ 



(29) 



and, 



Tu{z) = 4 J J d 3 r 1 d 3 r 2 G(f 1 - r 2 ) v(ri)v(r 2 )u(fi - z)u(f 2 - z)f\ ■ r 2 



+v(r 1 )v(f 2 - z)u(ri)v(f 2 - z)f 1 



[r 2 - z i 



\r 2 - z\ 



(30) 
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where u and v are, again, the upper and lower components of the single quark 
wavefunctions as in eq. 5. 

The spin sum, including the projection onto quark isospin, is given by 

g a P' a * = — 2 — ' ' 

so that the three body correction is given by 

5u 3b (x) = 25q 3b {x)-A8q 3b {x) 

5d 3b (x) = 6<$(x) + 6<&(x), (32) 



with 

x*)<,(*h + <M£)/ 3l, (-*), (33) 



and 

5 36f v _ £^ r d 3 k Mx-u-z-k 

xVW^h+VU^fM), (34) 

and where f 3b (k), f 3b {k) are the Fourier transforms of EB(z)f 3b (z) and 
EB(z)T 3b (z), respectively. Again the relative weights appearing in the eq. 
32 reflect the fact that a ud spectator pair is most likely to be in a spin zero 
state. 
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3 Gauge Correction 

As discussed in the last section, a gauge and renormalization scheme must be 
assumed in order for quark model wavefunctions to be meaningfully defined. 
Correspondingly, the effect of the path-ordered exponential of eq. 1 must 
be included in order to obtain quark distributions that are gauge invariant. 
Expanding the path ordered exponential P±(£) to leading order in g s , we 
obtain 

5q gau9 e(x) = ^jdcf o dr ] -(N ,P = dmC)A + (r ] -)^(0) 

x J d 3 rf(f) ■ A(f)\N , P = 0). (35) 

Neglecting self interactions, this expression may be evaluated in the same 
fashion as the two and three body corrections. The result, after some ma- 
nipulation, is given by 

x*P a (yP - g> + ^((i - y)P + q)z ■ Z p {yP - q, (1 - y)P + q). 

(36) 

The spin structure of the gauge correction is identical to that of the two-body 
correction, so we may immediately write 

Su gauge (x) = j-{2r^{x) + 2T?r 9e {x)) 

5d„(x) = ±(f^( x )+4Tr 9e (x)), (37) 
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where 

j gauge ^ 



16A 2 f 1 , fd 3 Pd 3 q^^ A . 

— tt \ dy / — — — —d (Mx — u) — z ■ q) 
fiM Jo J (2tt) 5 v h ' 

x Z l {k 2 , ki)z ■ k 2 (z ■ k 2 C w (ki, k 2 ) + z ■ k 1 C i(k 1 , k 2 ) 

+ (C 00 (k 1 , k 2 ) + h ■ k^Cnih, k 2 )) 
+Z 2 (k 2 , ki)z ■ k x {z ■ k 2 C w (k±, k 2 ) + z ■ k\C Q \{k^ k 2 ) 

+(C Q0 (k 1 , k 2 ) + h ■ hCuih, fc 2 )))l , 



(38) 



16A 2 



x 



3fiM Jo 



dy 



d 3 Pd 3 c 



5(Mx -u-z-q) 



x 



Z 3 (k 2 , ki)((ki ■ k 2 - z ■ hz ■ k 2 )C Q i(ki, k 2 ) 
-(k\ - (z ■ k 2 ) 2 )Cio(ki, k 2 ) + (z ■ k 2 k x ■ k 2 - k\z ■ k\)Cu{k u k 2 )) 
-Z 4 (k 2 , fci)((fci • k 2 - z ■ hz ■ k 2 )Cio{ki, k 2 ) 
-{k\ - (z • fci) 2 )C i(A;i, k 2 ) - {k\z -h-h- k 2 z • fci)Cn(A;i, k 2 )) 



(39) 



with ki = (1 — y)P + q and k 2 = yP — q. 

The gauge factor in eq. 1 removes the unphysical, gauge dependent phase 
in the quark wavefunction, which produces a mismatch between the momen- 
tum fraction x calculated using simple Fourier transforms and the physical 
momentum carried by the struck quark. This shift in x is manifested by the 
derivative with respect to x appearing in the gauge correction. This pat- 
tern persists in higher orders, as the path ordered exponential is expanded 
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in powers of g s -j~, and suggests that the gauge correction may be partially 



resummed to obtain a quark distribution of the form 



q\x) = (1 + 



d,Si(x) 



)q l (x + Si(x)), 



(40) 



dx 



where qb(x) is the zeroeth order quark distribution for quarks of flavor i. To 
leading order in a s , the shift function Si(x) is given by 



with Oj and bi the isospin dependent coefficients appearing in eq. 37. Since 
both a s and the derivatives of quark distributions can be large in the context 
of quark models, the effects of the resummation can be important, particu- 
larly for the d distribution. 

4 Results 

The expressions for the gauge, two and three-body corrections to the va- 
lence distributions of the proton were evaluated using a Gaussian quadrature 
scheme to perform the multiple integrations appearing in eqs. 20,21,33,34,38, 
and 39. For the gauge and two-body corrections, the function Zp{ki, k 2 ) was 
evaluated semi-analytically by making use of the 6 gaussian fit to the LAMP 
coordinate space wavefunctions described in reference 9. After accounting 
for the rotational symmetries of the problem, there were 4,5, and 6 integra- 
tions to be performed numerically for the two-body, gauge and three-body 
corrections, respectively. Treating all color magnetic effects perturbatively, 



Si(x) 



(q,/"(x) + ^r ge (x)) 
qi(x) 



(41) 
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we obtain the curves shown in Figs. 2 and 3 for the d and u valence distribu- 
tions, respectively. Also shown are the distributions with no color magnetic 
corrections applied, and the curves obtained by either neglecting the gauge 
correction, or by resumming the gauge correction in the fashion described at 
the end of section 3. 

As shown in Figs. 2 and 3, the general effect of CMI is to lower the 
momentum carried by the valence quarks. Physically, this makes sense, since 
the exchanged gluons will carry p + > 0, and consequently the momentum 
fraction carried by the struck quark must decrease. As advertised in the 
derivation of the two and three-body corrections, the d quark distribution is 
more sensitive to CMI than the u distribution, since the uu spectator pair is 
more likely to be found in a spin 1 state than the ud pair. Less encouraging 
is the fact that the Peierls-Yoccoz projection procedure is significantly less 
effective in producing corrected distributions with good support properties 
than it is for the unperturbed distributions. In particular, the two-body 
correction to the valence distributions has a non-negligible tail(~ .01) in the 
region x ~ 1. This increase in the relative size of the tail may be understood 
in terms of a small admixture of higher mass eigenstates in the Peierl-Yoccoz 
projected nucleon wavefunction. Since these states are more massive, the 
quarks they contain are not constrained to carry momentum p + < M / \[2 in 
the nucleon rest frame, and show up as a tail in the region x > 1. Since there 
are many excited states to couple to, the size of this unwanted component 
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of the wavefunction will tend to grow when a perturbation, such as the 



color magnetic interaction, is introduced JT2| . A more serious problem, from 
a phenomenological point of view, is the fact that the gauge correction is 
comparable to the two and three-body corrections. This suggests that the 
d quark distribution will be quite sensitive to the choice of gauge used to 
define quark model wavefunctions. In principle, this provides an additional 
restriction on the models, but in practice the infinite variety of gauges may 
provide a means to reconcile virtually any model with the data. 

The quark distributions are evolved from /i 2 , to 15 GeV 2 using a finite 
element procedure to reconstruct the evolved distributions from their mo- 



ments. A detailed description of this procedure will be given elsewhere fI5 
The results for the d and u distributions are shown in Figs. 4 and 5, respec- 
tively, for /Xq ranging from .2-. 4 GeV 2 . Also shown are BEBC neutrino data 
at the same Q 2 |I3 1 . In general, the data at large x favor a choice of the scale 
/zjj between .2 and .3 GeV 2 , while at smaller x, the data lies below model 
predictions even at these small renormalization scales. This situation may 
reflect the uncertainty associated with modeling the sea distributions, which 
are used to analyze the data in this region, or it may reflect a shortcoming of 
the model, namely that too much momentum is carried by the valence quarks 
at the quark model scale. In the former case, we note that parametrizations 



of parton distributions based on larger data setsJT^] generate larger valence 
distributions in the region of x in question. In the latter case, we have al- 
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ready noted that the momentum fraction carried by the valence quarks may 
be altered significantly by changing the spatial dependence of the empty bag 
matrix element, EB(z). 

In light of these facts, we conclude that color magnetic interactions pro- 
vide a natural mechanism for producing the observed differences between 
the u and d valence quark distributions of the proton, and that the LAMP 
picture provides a reasonable quantitative description of the valence distri- 
butions once color magnetic effects are included. 
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Figure Captions 



Figure 1 - Valence quark distributions for the MIT bag, soliton bag, 
and LAMP at the quark model scale. 

Figure 2 - Valence distributions for the d quark at the quark model 
scale, including color magnetic corrections. 

Figure 3 - Valence distributions for the u quark at the quark model 
scale, including color magnetic corrections. 

Figure 4 - Valence distributions for the d quark, evolved to Q 2 =15 
GeV 2 . Also shown are data from BEBC. 

Figure 5 - Valence distributions for the u quark, evolved to Q 2 =15 
GeV 2 . Also shown are data from BEBC. 
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